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Abstract 

Solar model independent tests of the vacuum oscillation and MSW solu- 
tions of the solar neutrino problem are considered. Detailed predictions for 
time (seasonal) variations of the signals in the future solar neutrino detec- 
tors (SNO, Super Kamiokande, BOREXINO, HELLAZ), if solar neutrinos 
take part in vacuum oscillations, are given. Results on the distortions of 
the spectra of neutrinos, and of e~ from the reaction u -\- e~ ^ u -\- e~ 
induced by ^B neutrinos, in the cases of vacuum oscillations or MSW transi- 
tions are presented for a large number of values of the relevant parameters. 
The possibilities to distinguish between the vacuum oscillation, the MSW 
adiabatic, and the MSW nonadiabatic transitions (solutions) in the future 
solar neutrino experiments are discussed. 
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1. INTRODUCTION 



With the pubhcation in 1991 and 1992 of the results of the Ga-Ge solar neutrino experi- 
ments [1,2] it became clear that the data from the first generation of solar neutrino detectors 
will not be sufficient to resolve the solar neutrino problem [3-6] which has been with us for 
more than 20 years. If the latest data provided by the pioneer Davis et al. [3] , Kamiokande 
[7], SAGE [1] and GALLEX [2] experiments are correct, an astrophysical explanation of the 
solar neutrino deficit seems unlikely at present [8] (for a recent discussion see [9] and [10]). 
At the same time, the current solar neutrino observations admit several rather different neu- 
trino physics interpretations which require the existence of unconventional neutrino intrinsic 
properties (mass, mixing, magnetic moment) and/or couplings (e.g., flavour changing neu- 
tral current (FCNC) interactions). These include: i) oscillations in vacuum [11] of the solar 
Vf. into different weak eigenstate neutrinos (z/^ and/or z/^, and/or sterile neutrinos, on 
the way from the surface of the Sun to the Earth [12] Q, ii) MSW transitions [13] Vf. — * '^/^{t); 
and/or Vf. Ug, while the solar neutrinos propagate from the central part to the surface 
of the Sun [14] 0, iii) solar Ue resonant spin or spin-flavour precession (RSFP) [16] in the 
magnetic fleld of the Sun [17], and iv) matter-enhanced transitions, for instance Ue Vr-, 
in the Sun, induced by FCNC interactions of the solar v^. with the particles forming the 
solar matter [18,19] (these transitions can take place even in the case of absence of lepton 
mixing in vacuum and massless neutrinos [18]) 0. Although the experiments of Davis et al., 
Kamiokande, SAGE and GALLEX will continue to run for at least few more years and the 

^For earlier discussions see, e.g., the references quoted in [12]. 

^The MSW solution has been studied by many authors before and after the publication of the 
results of the two Ga-Ge experiments: see, e.g., refs. [2,15] and the references quoted in [2,14,15]. 

^The solar neutrino decay hypothesis [20] is disfavoured [21], while mechanisms leading to universal 
suppression of the fluxes of ^B, ^Be, pp, etc. neutrinos due to Ve — > z^s transitions are ruled out, 
by the current solar neutrino data, if one uses the standard solar model (SSM) predictions of refs. 
[5,6] in the relevant analyses. 
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accuracy of the data they provided will improve, no substantial changes of the latter are 
expected Q and no qualitatively new data will be available before solar neutrino detectors 
of the second generation - SNO [22], Super Kamiokande [23], BOREXINO [24], ICARUS 
[25], and HELLAZ [26], become operational in the second half of the 90-ies 0. Our hopes 
for finding the cause of the solar neutrino deficit and for getting more precise information 
about the physical conditions in the central part of the Sun, where the neutrinos are being 
produced, are now associated with these future experiments. 

In the present article we continue the studies [14,27] of the possible solar model indepen- 
dent tests of the vacuum oscillation and the MSW solutions of the solar neutrino problem. 
The importance of these tests is difficult to overestimate given the fact that the solar model 
predictions for the neutrino flux may have rather large uncertainties. We present results 
on the specific seasonal time variations of the signals in the future solar neutrino experiments 
(SNO, Super Kamiokande, BOREXINO, HELLAZ), predicted if the solar neutrino deficit is 
caused by vacuum oscillations of solar neutrinos. We give also detailed predictions for the 
distortions of the spectra of *B neutrinos, and of e~ from the reaction u + e^ —* u + e" caused 
by the ^B neutrinos, in the cases of the vacuum oscillation and the MSW solutions. Neither 
seasonal time variations (apart from the standard ~ 7% geometrical one), nor substantial 
distortions of the spectra of the ^B, pp and CNO neutrinos (greater than ~ 10"^ E, E being 
the neutrino energy) are predicted to arise due to the specific physical conditions in the 
interior of the Sun [28]. The possibilities to distinguish between the vacuum oscillation, the 

priori, one cannot totally rule out the possibility of surprises in the next few years. The planned 
calibrations of the GALLEX and SAGE detectors will be crucial for the conclusive determination 
of the characteristics of the solar neutrino flux inferred from the current data. 
^Two of these detectors - SNO and Super Kamiokande, are under construction, BOREXINO and 
ICARUS are at the stage of prototype construction and/or testing, and the possibility to build 
HELLAZ is being studied. 
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MSW adiabatic, and the MSW nonadiabatic solutions of the solar neutrino problem using 
the data from the future SNO and Super Kamiokande experiments are discussed 0. Updated 
results on the MSW solution of the solar neutrino problem are also given. 

2. VACUUM OSCILLATIONS OF SOLAR NEUTRINOS: PREDICTED 
SEASONAL VARIATION EFFECTS AND SPECTRA DISTORTIONS 

The two-neutrino vacuum oscillation solution of the solar neutrino problem has been re- 
analyzed recently using the latest data from all currently operating neutrino experiments 



(Homestake, Kamiokande III, GALLEX and SAGE). The analysis was based on the predic- 
tions of the solar model of Bahcall and Pinsonneault ||^. It was found that the two-neutrino 
oscillations involving the Ue and an active neutrino, z/g ^ ^fiir), provide a not very good (but 
acceptable) quality of the x^^fit to the mean event rate solar neutrino data, while the oscilla- 
tions into sterile neutrino Ug, ^ I's, give a poor fit of the data: the u^. ^ '^^i{t) oscillations 
are ruled out at 90 % C.L., but are acceptable at 95 % C.L., while the z/g z/g oscillations 
are ruled out at 99 % C.L. The results are rather different if one uses the data available 
from each particular run of measurement of the Homestake, Kamiokande II, GALLEX and 
SAGE collaboration in the x^^analysis (for details see ref. ||30|). Both the z/g ^ z/^(^) and 



^ '^s oscillation hypotheses give good fits to the run-by-run solar neutrino data, being 
acceptable even at 68 % C.L. The regions of values of the two parameters, Am^ and sin^ 29, 
characterizing the two-neutrino oscillations of the solar z/g, which are allowed (at 90 % C.L.) 



by the run-by-run data, lie in the following narrow intervals [pO| : 

i^e ^ z^mW : 5.7 X 10""eV^ ~ Am^ ~ 1.1 x 10~%V^ (la) 

sin2 20~O.75, (16) 



^For alternative solar model independent tests of these solutions, based on the future SNO and 
Super Kamiokande data see refs. [29]. 
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and 

z/e ^ z/, : 5.0 X lO^^eV^ ~ Am^ ~ 6.6 x IQ-^eV^ (2a) 

sin2 2^~0.8. (26) 

2.1 Seasonal Variations of Signals 

The probability that a solar electron neutrino with energy E will not change into z/^(t-) 
(or Us) on its way to the Earth when z/g z/^(t-) (z/g ^ z/^) oscillations take place, has the 
form: 

P(z/e ^ z/e;R(t),E) = 1 - -sin2 2^ [1 - cos27r5^], (3) 

2 Lv 

where Ly = 47rE/Am^ is the oscillation length in vacuum, 

R(t) = Ro [l-ecos27r;i], (4) 

is the Sun-Earth distance at time t of the year (T = 365 days), Rq = 1.496 x 10^ km and 
e = 0.0167 being the mean Sun-Earth distance and the ellipticity of the Earth orbit around 
the Sun. 

For E ^ 1 MeV and the values of Am^ from the intervals (la) and (lb) one has: 
Lv — 27r(2eRo), where 2eRo is the variation of the Sun-Earth distance in the period De- 
cember - June. This implies that if solar neutrinos take part in vacuum oscillations, the flux 
of solar neutrinos will exhibit seasonal variations []. The magnitude of the time variations 
depends, in particular, on the energy of solar neutrinos and will be different for the ^B, 
^Be, pp, pep and the CNO neutrinos. Obviously, if the integration over the neutrino energy 
renders the oscillating term in the expression for P(z/e z/e;R(t),E) negligible (as is the 



case of pp neutrinos [0]), the energy integrated observables will not exhibit seasonal (time) 
variations. 



^The possibility of seasonal variations of the flux of solar neutrinos when the latter take part in 
oscillations in vacuum was indicated first by B.Ya. Pomeranchuk (see, e.g., ref. ||3l[| ). 
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We have depicted in Figs. 1 and 2 the expected time variation of the ratio of the predicted 
signals (event rates) due to the ^B, pp, ^Be and pep neutrinos in the case of (two-neutrino) 
oscillations of solar neutrinos, to the ones in the absence of oscillations. The results shown 
in these figures are valid for experiments detecting the solar neutrinos via the u — e~ elastic 
scattering reaction (Super Kamiokande, BOREXINO, HELLAZ, etc.). The electron kinetic 
energy detection threshold was taken in the calculations of the signals due to the ^B, the 
pp, the ^Be and pep neutrinos to be 5 MeV, 0.1 MeV, and MeV, respectively; the possible 
effects of the detectors efficiencies and finite electron energy resolution were not included in 
the calculations (we leave it to our colleagues-experimentalists to take into account these 
effects in accordance with the specific characteristics of their respective detectors) . Electrons 
with kinetic energy Tg > 5 MeV will be detected in the Super Kamiokande experiment. The 
HELLAZ detector is planned to be sensitive to e~ (from the pp neutrino induced reaction) 
with Te > 0.1 MeV. As for the BOREXINO detector in which ~ 90% of the event rate 
is predicted to be generated by the 0.862 MeV ^Be neutrinos, it is expected that the signal 
to background ratio will allow one to extract the ^Be neutrino signal for e~ with kinetic 
energy in the interval 0.25 MeV < Tg < 0.66 MeV. We have checked that reducing the 
interval < Tg < 0.66 MeV used in our calculations to the one refiecting the currently 
envisaged detection capabilities of BOREXINO has no observable effect on the results (the 
two ^Be(862) curves, corresponding to the two different intervals of integration would be 
indistinguishable if both were plotted in Figs. 1 and 2). 

In the calculations of the signals (event rates) in the case of vacuum oscillations we 
have taken into account also the standard R~^(t) dependence of the values of the different 
solar neutrino fiux components at the Earth surface. Finally, the SSM predicted event 
rates (signals) have been obtained by dividing by T the one year total number of events 
calculated within the SSM [5] (assuming 100% detection efficiency). Thus, in the absence of 
vacuum oscillations the plotted (theoretical) ratios will change from (1 — e)~^ in December 
to (1 + e)~^ in June, while the ratio of measured to the SSM predicted [5] signals will vary 
from A(l — e)~^ to A(l + e)~^, where A is a constant which can be different for the different 
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(pp, ^Be, ^B, and pep neutrino induced) signals (A = 1 if the SSM prediction [5] for the flux 
of the corresponding neutrinos (pp, and/or ^Be, and/or ^B, and/or pep) is correct). 

The results shown in Figs. 1 and Figs. 2a-2d correspond to solar z/g oscillations into 
active neutrino, Vf, <-> z/^(t-) , while in Fig. 2e and Fig. 2f we have depicted results in the case 
of oscillations into sterile neutrino, v^. 

As was indicated on the basis of few numerical examples in [14,27], the most dramatic 
seasonal variations are predicted to be exhibited by the signals due to the monochromatic 
^Be and pep neutrinos. Typically, the differences between the signals in December and June 
are the largest. However, as was noted in for certain values of the parameter Am^ the 
signals in December (or June) and March (or September) differ most (see, e.g.. Figs. Ic, Id 
and 2f). Figs. 1 and 2 demonstrate also that the predicted magnitude and explicit form of 
the time variations of the ^Be and pep neutrino induced event rates in the v — e~ elastic 
scattering experiments are extremely sensitive to the value of Am^ and change drastically 
even for relatively small variations of this parameter. 

Since the pp neutrinos have a rather low energy (E < 0.42 MeV), for most of the values of 
Am^ from the intervals (la) and (2a) the inequality l-nRo » holds. As a consequence, 
the integration over the recoil e~ energy in the calculations of the pp neutrino induced signals 
renders the oscillation term in the probability (3) negligible [|l^. Therefore the predicted 
seasonal change of the energy integrated signals due to the pp neutrinos in the u — e~ 
scattering experiments coincides with the standard 7% geometrical one (Figs. 1, 2a-2c, and 
2e), except for values of Am^ ^ 5 x 10^^^ eV^, for which the vacuum oscillations lead to 
rather small deviations from it (Figs. 2c and 2f). 

The seasonal changes of the signal due to the ^B neutrinos do not exceed approximately 



15% Such variations are not detectable in the currently running experiments. However, 
the high statistics future solar neutrino experiments Super Kamiokande, SNO, and ICARUS 
are envisaged to accumulate (between 3000 and 4000 events per year) will allow to detect 
even rather small (few percent) differences between the signals in December and June. For 
Super Kamiokande (ICARUS) and SNO detectors the effect of time variations is shown 



separately in Figs. 3a (3b), 3d (3e), and in Figs. 3c and 3f, respectively, using a different 
normalization of the signals and a proper scale. The predictions for the Super Kamiokande 
detector depicted in Figs. 3a (3d) and 3b (3e) differ in the value of the neutrino threshold 
energy, Eth, used in the calculations (see further): Eth — 5 MeV and Eth — 7.5 MeV, 
respectively. The smaller magnitude of the effect in comparison with that in the case of the 
signals, generated by the monoenergetic ''Be and pep neutrinos, is not difficult to understand 
quahtatively. The Super Kamiokande and SNO experiments will be sensitive only to 
neutrinos having relatively high energies (E > 5 MeV and E > 6.44 MeV) which exceed at 
least by a factor of 6 the energy of the dominant (0.862 MeV) component of the ''Be neutrino 
flux. For these energies and the values of Am^ from the intervals (la) and (lb) one has: 
27r(eRo/L^) < 0.14. As it follows from eqs. (3) and (4), under this condition the seasonal 
changes of the probability P{i'e ~^ ^^e;R'(t))E) are proportional to, and do not exceed, the 
ratio 27r(eRo/L^) and, therefore, cannot be large. The integration over the neutrino energy 
reduces further the magnitude of the effect of interest. 

The solid, dotted, dashed, long-dashed, dash-dotted, and long-dash-dotted lines in 
Figs. 3a, 3b and 3c (Figs. 3d, 3e and 3f) represent results for the values of Am^ and 
sin^2^ for which Figs, la, lb,.... If (Figs. 2a, 2b,..., 2f) have been obtained. The nor- 
malization of the signals shown graphically in Figs. 3 is chosen in such a way as to avoid 
any dependence on the prediction for the total flux of ^B neutrinos, and thus on the so- 
lar models. Namely, for given Am^ and sin^ 26 the calculated event rate at time t of the 
year in the case of vacuum oscillations, dNe^{Am'^,9,Eth',t)/dt, is divided by the quan- 
tity dNl{t)/dt = Ne„(Am2, e, Eth, ly) T-i(Ro/R(t))2, where Ne^(Am2, 6, Eth, ly) is the pre- 
dicted total number of events per year provided ^B neutrinos undergo vacuum oscillations 
with the chosen values of the parameters Am^ and sin^ 26, and the ratio 

2 dNe,(Am^g,E,ft;t)/rft R^{t) dNe,(Am^ E^^; t)M 

dKviV/dt Ro Nev{Am^,6,Eth;ly) 
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is plotted in Figs. 3. For the SNO detector we have: 



dNev{Am^,e,Eth;t) 1 



14.4 MeV 



dt R2(t) 



J FSS^n(E) P(i/e ^ i^e; R(t), E) (7(i/ed ^ e-pp) 
Ne^(Am , 9, Eth] ^y)^ J dt, (7) 



where F^^'^/R^(t) is the predicted total flux of ®B neutrinos at the Earth surface at time 

14.4 MeV 

t of the year, n(E) is the normalized to 1 spectrum of neutrinos, / n(E)d£' = 1, 



Etfi = 6.44 MeV, and a{ved — ^ e^pp) is the cross-section of the charged current reaction 
i/g + d ^ e~ + p + p hy which the solar neutrinos will be detected in the SNO experiment. 
Obviously, expression (6) corresponds to ideal detection conditions; for the comparison of 
the theoretical predictions with the future SNO data it has to be modified by taking into 
account the neutrino energy resolution function, the detection efficiency, etc. of the SNO 
detector. The expression for the predicted event rate in the Super Kamiokande detector can 
be obtained from eq. (6) by replacing alu^d —>■ e~pp) with the cross-section cr(fee~ — > t'eC") 
of the reaction i/g + e~ — > + e~, and by using an appropriate value for Eth] in the case of 
^ij.{t) oscillations the probability P(fe — > ^e] R-lt), E) must be substituted with 

r, + (l-r,) P(i/e^i^e;R(t),E), (8) 

where — (T(^'^e~ v^e~^la{ye.e' — > v^e' 



~ 1 

~ 6" 



It is not difficult to convince oneself that c?N(t)g^/di is the event rate at time t of the 
year if the total number of events per year is Ne^(Am^, E^/^; 1 y) and the relevant signal 
does not exhibit any additional time dependence, except for the standard R~^(t) geometrical 
one. The ratio (5), evidently, is independent of the value of F^^'^ and thus is solar model 
independent. The comparison of the predictions presented graphically in Figs. 3 with the 
data will be straightforward: as input one needs only the experimentally measured mean 
event rate for a given interval of time (one month, say), and the total number of events 
observed per year; the latter will provide the value of Net,(Am^, Et/^; 1 y). All the other 
quantities entering into the ratio (5), T, Rq and R(t), are known with a high precision. In 



the absence of vacuum oscillations the ratios of signals (5) plotted in Figs. 3 will be equal 
to 1. 

We shall obtain next an approximate but sufficiently accurate and rather simple analytic 
expression for the time variation observable (5), exhibiting its time dependence explicitly. For 
e = 0.0167 < 1 and 27r(eRo/L„) < 0.14 < 1, the quantity 27r(eRo/L„) cos(27rt/T) entering 
into the formula for the probability P{i'e i^e; R-lt), E), can be used as a small expansion 
parameter together with ecos(27r;|;). Expressing the oscillating term in P(fe — > fe;R(t),E) 
as a power series in 27r(eRo/L^) cos(27rt/T), and R(t) as a power series in ecos(27r7|), it is 
easy to show that the leading correction in Ne^(Am^, 9, Eth', 1 y) due to the ellipticity e is 
proportional to and does not exceed 5 x 10~^. Thus, up to corrections ~ 5 x 10~^, the 
quantity Ne^(Am^, 9, Efh] 1 y) does not depend on e and can be obtained by setting e to in 
eqs. (6) and (7). Using this fact, and keeping in (6) only the terms up to the second order in 
27r(eRo/L^) cos(27rt/T) in the expansion of P{ue — > Ue', R(t), E), one arrives at the following 
result for the observable Rvar(Am^, 9, Eth;t) for the SNO detector: 

R^O(Am^e,Eth;t) = l + ecos(27r;^) sin'29 K^''''{Am',9,E,^), (9) 

where 

14.4 MeV 

/ dE x(sin 2x — xcos 2x ecos 2tt^) n(E) a{i/ed — > e~pp) 

K'^^{Am',9,E,^) = +0{{2xn 

J dE P(i/e Ve] Ro, E) n(E) a{v^d e-pp) 

(10) 

and X = evrRo/Lv < 0.07. The corresponding expression for the time variation observable 
for the Super Kamiokande detector, R^(Am^, 9, Eth; t), can be obtained formally from eqs. 
(9) and (10) by replacing in eq. (10) the probability P(i/e — i^e;R(t),E) by [rj, + (1 — 
ij,) P(z/e ^e] R(t), E)], al^ed —>■ e~pp) with (T(z/ee~ — > t'eC"), by changing the value of Eth, 
and by multiplying the numerator in eq. (10) by the factor (1 — r,^). Let us note that the 
effect of the time dependence of KSno(sk)(^^2^ g^^^ R^o(sk)(^^2^ g^^. -g beyond 

the sensitivity of the next generation of experiments. 

Few comments concerning the results shown in Figs. 3 are in order. All (correspondingly 
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normalized) signals are equal to 1 at t = |T and t = |T, in accordance with eqs. (9) and 
(10). As Figs. 3 indicate, the predicted amplitude of the time variations of the signal in the 
SNO detector is typically (but not always, e.g., compare the dashed and the long-dashed 
hues in Figs. 3a and 3c) larger than that in the Super Kamiokande detector. The difference 
in the magnitude of the signal time variations in the two detectors is a consequence of i) the 
difference in the minimal ®B neutrino energy the two detectors are planned to be sensitive 
to (6.44 MeV and 5 MeV), ii) the specific neutrino energy dependence of P(i/e — ^ T^e] R'(t), E) 
in December and June, iii) the difference in the El-dependence of the cross-sections a^u^d 
e~pp) and a{i'ge~ — > t'eC"), and in the case of oscillations iv) the probability 

P{ue — > Ue', R'(t), E) entering into the expression for the predicted signals in SNO and Super 
Kamiokande detectors with different coefficients (see eqs. (6) and (8)): 1 and approximately 
|, respectively. These differences can lead even to a strong anticorrelation between the 
signals in SNO and Super Kamiokande experiments, as in the case of Am^ = 6.9 x 10~^^ eV^ 
and sin^ 2^ = 0.9 (the sohd lines in Figs. 3d and 3f). 

If |K^'^°(^^)(Am^, 6', Eth)| -C 1 for certain values of Am^ and sin^ 2^, one has 
R|m°^^^^ (Am^, 9, Eth; t) = 1+0(10"^), and the time variation effects will not be observable in 
SNO (Super Kamiokande) experiment in spite of vacuum oscillations of neutrinos. Such 
is practically the case with the signal in the Super Kamiokande detector for E^h = 5 MeV, 
Am^ = 6.3 X 10-1^ eV^ and sin 29 — 0.85 (see Figs, le and 2e and the dash-dotted lines 
in Figs. 3a and 3d). Fortunately, our results show that the indicated possibility is never 
realized both for the signals in the Super Kamiokande and the SNO detectors (compare the 
dash-dotted lines in Figs. 3a, 3d and in Figs. 3c and 3e). Moreover, it can take place either 
for the event rate in the Super Kamiokande detector measured with Eth — 5 MeV, or for 
the event rate obtained with = 7.5 MeV (the long-dashed line in Fig. 3b), but not for 
both event rates (compare the dash-dotted lines in Figs. 3a, 3d and in Figs. 3b, 3e, and the 
long-dashed lines in Figs. 3b and 3a). 

In certain cases the magnitude and the pattern of the time variation of the signal in 
the Super Kamiokande detector is very sensitive to the increase of the threshold neutrino 
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energy from Eth = 5 MeV to Eth = (7 — 8) MeV. This is illustrated in Figs. 3b and 3e, 
where the results of the calculations of the ratio of signals (5) for the same values of Am^ 
and sin^ 29, for which Figs. 3a and 3d were obtained, but with Eth = 7.5 MeV (instead of 
Eth = 5 MeV), are presented. We see, in particular, that for Am^ = 6.9 x 10^^^ eV^ and 
sin^ 26 = 0.9 the pattern of the time variations has changed completely with the change of 
Eth'- the maximum of the ratio (5) is now in December rather than in June, and the predicted 
variations of the signals in the SNO and Super Kamiokande detectors are correlated (rather 
than anticorrelated). Depending on the value of Am^ (and sin^2^), the change of Eth from 
5 MeV to 7.5 MeV can increase, or diminish the amplitude of the variations (compare, e.g., 
the solid, dotted, dash-dotted, and long-dash-dotted lines in Figs. 3a and 3b, as well as the 
dashed and long-dashed lines in the same two figures); for some values of Am^ and sin^ 20 
the increase is quite substantial. We have not studied the effect of change of Eth on the time 
variation of the signal in the SNO detector. However, one can expect on the basis of the 
above results that for certain values of Am^ and sin^ 29 it can be dramatic. 

It is interesting to note also [14,27] that for certain values of the parameters Am^ and 
sin^ 29 the seasonal change of the neutrino induced signals, associated with the vacuum 
oscillations, can compensate partially or completely the standard 7% geometrical one and 
in the second case the event rate (iNev(Am^, 0, Eth; t)/(it will be constant in time (see Figs. 
la-Id and 2a and the corresponding solid, dotted, dashed and long-dashed lines in Fig. 3a, 
the dotted line in Fig. 3b, the dashed and long-dashed lines in Fig. 3c, as well as the solid 
line in Fig. 3d); it can even lead to an increase of the event rate dNey{Am'^,9,Eth',t)/dt 
from December to June (see the solid line in Fig. 3b and the solid and dotted lines in 
Fig. 3c). Note that due to the specific normalization chosen by us a constant event rate 

^In the case of the monoenergetic ^Be and pep neutrinos even a dramatic increase of the corre- 
sponding signals from December to June due to the vacuum oscillations is possible (as can be seen 
in Figs. 1 and 2). 



12 



will correspond to an increase of the ratio (5) plotted in Figs. 3 from the value (1 — 2e) in 
December to the value (1 + 2e) in June; an increase of the event rate from December to June 
corresponds to an increase of the ratio (5) from a value smaller than (1 — 2e) in December to 
a value greater than (1 + 2e) in June (see the solid line in Fig. 3b and the solid and dotted 
lines in Figs. 3c). Thus, a non-observation of the 7% change of the neutrino induced 
event rate (constant rate), or a registration of an increase of the rate, in the period from 
December to June in SNO and/or Super Kamiokande detector would be a strong evidence 
that solar neutrinos take part in vacuum oscillations. Note that, as is clear from Figs. 3, for 
given [S.W? and sin^ 2Q the compensation (partial or complete) of the standard 7% seasonal 
variation can take place either for SNO or for Super Kamiokande signals, but not for the 
signals in both detectors. Futhermore, in the case of the signal in the Super Kamiokande 
detector such a compensation does not hold both for E^/i = 5 MeV and for Ej/j = 7.5 MeV. 

2.2 Spectra Deformations 

If solar neutrinos take part in vacuum oscillations, the shapes of the spectra of the ^B, pp, 
and the CNO neutrino fluxes at the Earth surface will differ from their standard forms. The 
corresponding spectra deformations will reflect the specific and relatively strong dependence 
of the oscillation probability P(z/e ^'o;R(t),E), eq. (3), on the neutrino energy E. The 
change of the solar neutrino spectrum will lead also to a change in the spectrum of the final 
state electrons in the z/ — e~ elastic scattering reaction induced by the solar neutrinos. 

The deformation of the (average) spectrum of *B neutrinos for the same 12 values of 



The spectra under discussion will also exhibit relatively small seasonal variations if solar neutrinos 
undergo vacuum oscillations. Here we have in mind the average spectrum which will be determined 
experimentally from data collected during a period of k years, k = 1,2,3,... . The relative magnitude 
of the correction due to the seasonal variations in the average spectrum of neutrinos is not greater 
than ~ 5 X 10^^, while the relative difference between the spectra in December and June does not 
exceed 14%. 
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the parameters Am^ and sin^ 26, for which Figs, la-lf, 3a-3c and Figs. 2a-2f, 3d-3f have 
been obtained, are shown respectively in Fig. 4a and Fig. 4b. Each (average) spectrum, 
(i$B(Am^, ^, E)/(i£', to be determined from data collected by the SNO experiment over a 
period of k years, 

(i$B(Am^g, E) _ I 1 ciNev(Am^g,E;k) _ 

dE kT aiPgd — ^ e~pp) dE 

pSSM 

= — ^ n(E) P(i/e ^ ^/e;Ro,E), (11) 

Kg 

(iNev(Am^, 9, E; k) / dE being the total number of events induced by ®B neutrinos with energy 
E in k=l,2,3,... years, while the last term in eq. (11) represents the theoretical expression for 
the spectrum in the case of vacuum oscillations, is divided by the (average) SSM spectrum, 

t;$|SM(E) fssm 



dE Rl 



n(E), (12) 



predicted in the absence of oscillations. To avoid the dependence on the SSM prediction for 
the total flux of neutrinos this ratio of spectra, 

rSNO./\^2 n ^^ _ ^^B(Am^g,E)/cii? 

is further normalized to the value of the ratio at E = 10 MeV, and the double ratio 
RSpNO(Am2, e, E) (n(E) a{uA -> e-pp))E=io MeV d^,,{/\ui\ 6, E; k)/c/E 



RSpNO(Am2, e, 10 MeV) n(E) a{uA ^ e'pp) dNev(Am2, 9, 10 MeV; k)/rfE 

^ P(z/e ^ t^e; Rq, E) 

P(i/e ^ i^e;Ro,E= 10 MeV)' ^ ' 

is plotted in Figs. 4a and 4b. Thus, in the case of absence of deformations the ratio of 
spectra depicted will be constant (i.e., neutrino energy independent) and equal to 1 Note 



^'^The absolute deformations of the spectra of the and pp neutrinos in the case of i>e ^ ^ii{t) {^^ 
Ve <-> i^s) oscillations, and for the SSM predictions of ref. [5] have been shown in ref. [27] for four 
pairs of values of the parameters Am^ and 'sit?2Q, namely, for (Am2[eV2]; sin^ 20) = (1.1x10^^^; 
1.0), (7.9x10-11; 0.8), (6.3x10""; 0.8), (5.5x10""; 1.0). 
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that this would be vahd both for a constant reduction of the spectrum of the flux (and 
therefore of the total flux) of '^B neutrinos by a certain (energy independent) factor, and if 
there is no reduction at all and the flux coincides with the predicted one. 

The changes of the (average) spectrum of the final state e~ in the z/— e~ elastic scattering 
reaction induced by the neutrinos in the cases of z/g ^ ^^l{T) Q ^ oscillations 

are shown respectively in Fig. 5a, Fig. 5b (curves labelled 1-4), and in Fig. 5b (curves 
labelled 5 and 6), for the same 12 values of Am^ and sin^ 29 for which Figs. 1, 2, 3 and Figs. 
4a, 4b have been obtained. The e~ kinetic energy range chosen (5 MeV < Tg < 14 MeV) 
coincides with the one to which the Super Kamiokande detector is planned to be sensitive. 
The recoil-electron spectra depicted in Figs. 5a and 5b are normalized in the same way as 
the spectra shown in Figs. 4a and 4b Q, i.e., the following double ratio is plotted in Figs. 
5a and 5b: 

R3SpK(Am2, 6, Te) _ ^ dN,,{Am^, 9, T,; k)/dT, 



w 



(T 



RSpK(Am2, 6, 10 MeV) ' ' dN^,{Am^, 6, 10 MeV; k)/dTe 

14.4 MeV 

/ n(E) (il + (1 - T^Fiu, ^ z/e; Ro, E)) {daiu^e- ^ Vee-)/dT,) dE 

Tc(l+^) 

= ^(Te) 14.4 MeV ' (1^) 

/ n(E) (r', + (1 - r'jP(z/e ^ Ro, E)) {daiv^e- ^ Vee-)/dT,) dE 

10.25 MeV 

where (iNev(Am2, 9, Tq; \C)/dT(. is the number of events (observed in k years) with the recoil 
e~ having an energy Tg, da{v^e~ — > Ve^~)/dT^ is the differential cross-section of the process 
z/g + e" ^ z/g + e", = {da{y^e' Uf,e~) / dTe) / {da{uee~ Uee~)/dTe) = (| - y), and 



^^More precisely, induced by the "surviving" electron neutrinos and by the i^ij,(t) neutrinos into 
which the neutrinos have oscillated. 

^^In ref. (see Fig. 3a) we have shown just the ratio of the e~spectrum in the case of fg ^ ^^{t) 
oscillations, and of the standard e~spectrum, for four pairs of values of Am^ and sin^ 29, chosen 
from different parts of the intervals (la) and (lb). 
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14.4 MeV 

/ n(E) {da{vee- v^e-)ldT^ dE 

w(T.) = ^etv ■ ( 16) 

/ n(E) ((i(j(i/ee~ —>■ i'ee~)/dTg) dE 

Te(l+^) 

Thus, in the absence of deformations (no reduction, or energy independent reduction of the 
electron neutrino flux) the double ratio of e~— spectra (15) will represent a horizontal 
line crossing the vertical axis at the point 1. 

Let us note that one can choose to normalize the ratios of the predicted and the standard 
spectra discussed above by their values not at 10 MeV, but at some other (in general, different 
for SNO and Super Kamiokande detectors) energies. For a given experiment the energy of 
normalization must be chosen on the basis of considerations of accuracy of the corresponding 
data, and of maximal enhancement of the effect of deformation if present in the spectrum. 

One can utilize an alternative spectrum normalization based on the measurement of the 
total (average) flux of neutrinos with energy E> Eth to form a solar model independent 
observable. In the case of the SNO detector this total flux is given by 

14.4 McV r;^ / ^ 9 /I T^x 

rf<l>B(Am2,^,E) 

where the integrand is determined by eq. (11). In the absence of vacuum oscillations (or 
MSW transitions) the spectrum of ^B neutrinos having energies E> Eth, whose total flux 
is $B(Am2, e, Eth), will have the form: d^^{E)/dE = n(E) *B(Am2, 6, Eth). The total flux 
$B(Am^, 9, Eth) (or the spectrum d^^{E)/dE) can be used to normahze the measured spec- 
trum (11). Thus, instead of the double ratio (14) one can consider solar model independent 
ratio 

rSNO/ .^2 n _ 1 d$B(Am^g,E) _ 
(^"^'^'^)-$B(Am^^,Eth) dE 

_ (a(t/ed ^ e-pp))-^ dNev(Am^ 9, E;k)/dE _ 

~ 14.4 MeV 

i dE {a{v^d e-pp))-i dNev(Am2, 9, E; k)/d£; 

Eth 

= n-(An.'.g.E.,) P(--. ^ R^. E). (18) 
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where 

14.4 MeV 

n'{Am\ e, Eth) = J n(E) P{i^, ^ p,- Rq, E) dE (19) 

Eth 

is the total neutrino flux suppression factor in the case of vacuum oscillations, < 
n'(Am^, ^, Eth) < 1- The analogous ratio for the Super Kamiokande detector can be easily 
derived. There are two advantages in utilizing the normalization described above: i) the 
corresponding ratios of spectra will de determined experimentally with a higher precision 
than the double ratios (14) and (15), and ii) it allows a straightforward comparison between 
the theoretical predictions and the data. For certain values of Am^ and sin^ 29 the spectra 
deformations can be less pronounced in the ratios of the type (18) than in the double ratios 
(14) and (15), and vice versa. This is iUustratcd in Figs. 5a and 5b, where we show the 
®B neutrino spectra depicted respectively in Figs. 4a and 4b, but normalized in the manner 
described above, eq. (18). 

3. MSW TRANSITIONS: IMPRINTS ON THE SPECTRA 

In the case of two-neutrino MSW transitions — > or i/e — > J^s in the Sun, the 
solar Ve survival probability, P(t'e — > Ve] E), can be calculated with very high accuracy for 
Am^ ^ 5 X 10~^ eV^ and sin^ 29 > 10~^ using the simple analytic expression [32,33]: 

P{Ue ^Ue-, E) = ^ + {^- P') COs2^„(to) COS 2^. (5) 



Here 



^, exp[-7rro ^(l-cos2g)]- exp[-27rro-^™'i 
1 - exp[-27rro- 



is the level crossing probability (i.e., the analog of the Landau-Zener probabihty) for the 
case of density varying exponentially along the neutrino trajectory in the Sun, 9m{to) is the 
neutrino mixing angle in matter [13] in the point of Vg production in the Sun, and ro is the 
"running" scale height [32,33] (see also [14]), i.e., the scale height calculated at the resonance 
point. For Am? > 5 x 10^^ eV^ and sin^ 29 > 10~^ expression (5) allows one to derive MSW 
predictions for any observable quantity associated with the detection of the solar neutrino 
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flux, or of its different components, on Earth 0. 

We have re-examined (exploiting the x^— method) the MSW solution of the solar neutrino 
problem using the most recent published data from all four operating solar neutrino detectors 
(see Fig. 7). The analysis was based on the SSM predictions of ref. [5]. It revealed that 
in the case of z/g ^h{t) transitions i) the "lower" (in values of Am^) branch of the large 
mixing angle (adiabatic) solution [14] (actually, the region /S.rn? < xlO~^ eV^, sin^ 2^ > 0.1) 
is excluded by the current data at 99.5% C.L., ii) the "upper" branch [14] provides a not 
very good quality of the fit of the data (min = 5.30 (with the theoretical uncertainties 
included in the analysis) for 2 d.f.), being excluded at 90% C.L., but allowed at 95% C.L., 
and iii) the small mixing angle (nonadiabatic) solution [14] provides the best fit of the data 
(min = 0.48 (with the theoretical uncertainties included in the analysis) for 2 d.f.). The 
results are quite different if one assumes that z/g Ug transitions take place: in this case 
only a small mixing angle (nonadiabatic) solution is acceptable (at 90% C.L.: min = 3.43 
(without the inclusion of the theoretical uncertainties in the analysis) for 2 d.f.), while a 
large mixing angle solution is excluded at 99.7% CL. Note, however, that the nonadiabatic 
solution in the case of ^h(t) transitions gives a better quality of the fit of the data 

than the nonadiabatic solution associated with the z/g transitions. Our results are 

presented graphically in Figs. 7, where the regions of values of Am^ and sin^ 29, allowed at 
90% CL. and 95% CL. are depicted: Figs. 7a, 7c and Figs. 7b, 7d correspond respectively 
to z/g z/^(t-) and z/g — conversions. 

The distortions of the spectrum of the neutrinos (E> 5 MeV) predicted in the 



very precise and simple analytic description of the two-neutrino MSW transitions of solar 
neutrinos for sin^ 26* < 10^'^ was derived in ref. [33]. If Am^ ^ 5 x 10^*^ eV^, for sin^ 20 > 0.1 
one must take into account in the description of the transitions of the monoenergetic ^Be and pep 
neutrinos also the nonadiabatic oscillating term present in P(t'g Ve] E) [34], for which there 
exists a relatively simple analytic expression as well [32]. 
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case of two-neutrino MSW transitions are shown in Figs. 4c, 4d, and 5c, 5d, while the 
corresponding distortions of the spectrum of e~ from the reaction v -\- e~ ^ v -\- e~ induced 
by the neutrinos are depicted in Figs. 6c and 6d. The MSW spectra shown in Figs. 4c 
(5c), 4d (5d) and 6c, 6d are normalized in the same way as the vacuum oscillation spectra 
depicted in Figs. 4a (5a), 4b (5b) and 6a, 6b. Thus, plotted in Figs. 4c (5c), 4d (5d) and 6c, 
6d are the corresponding double ratios (14) and (15) (ratio (18)). In previous publications 
we have shown graphically just the ratio of the predicted MSW and the standard e~— spectra 
(ref. [14], Fig. 3b), and the absolute deformations of the and pp neutrino spectra (ref. 
[27], Fig. 5d) for the same four pairs of values of Am^ and sin^ 26*, for which Figs. 4c and 
6c are obtained. As is evident from Figs. 4-6, the measurements of the ^B neutrino and of 
the recoil-electron spectra in SNO, Super Kamiokande and ICARUS experiments will allow 
one, in particular, to discriminate between the MSW nonadiabatic and the MSW adiabatic 
solutions of the solar neutrino problem. 

4. DISTINGUISHING BETWEEN THE VACUUM OSCILLATION AND 
THE MSW SOLUTIONS 

An unambiguous evidence of vacuum oscillations of solar neutrinos would be the observa- 
tion of clear deviations from the standard 7% seasonal variation of the signals in the future 
solar neutrino detectors: no other solution of the solar neutrino problem leads to such an 
effect. In the case of vacuum oscillations the predicted nonstandard seasonal changes of the 
signals due to the monoenergetic ^Be and pep neutrinos are the most dramatic (see Figs. 1 
and 2). Although much smaller, the seasonal variation effects in the signals generated by the 
^B neutrinos are, for most of the values of Am^ and sin^ 29 from the intervals (la) and (lb), 
sufficiently large to be detected by the SNO, Super Kamiokande and ICARUS experiments, 
provided the detectors will operate with their envisaged detection capabilities and expected 
background levels. As we have demonstrated, the effects can be enhanced by choosing ap- 
propriate values of the relevant threshold detection energies. The data on the seasonal time 
variations of the event rates in SNO, Super Kamiokande and ICARUS experiments can be 
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crucial for discriminating between the vacuum oscillation and the other possible solutions of 
the solar neutrino problem. 

The predicted distortions of the *B neutrino and the recoil-electron spectra due to two- 
neutrino vacuum oscillations or MSW transitions of the neutrinos (Figs. 4, 5 and 6) 
provide us with an indispensible possibility to test these solutions in a solar model indepen- 
dent way in SNO, Super Kamiokande and ICARUS experiments. As is evident from the 
comparison of Figs. 4a,..., 4d (5a,..., 5d) and Figs. 6a,..., 6d, respectively, both the vacuum 
oscillations and the MSW transitions lead to somewhat stronger shape deformations of the 

neutrino spectrum than of the recoil-electron spectrum: some of the features of the dis- 
torted ^B neutrino spectrum are less pronounced, or are not present, in the e~ spectrum as a 
result of the integration over the neutrino energy necessary to perform to obtain the latter. 
The only exception are the spectra corresponding to large mixing angle MSW transitions 
(see the curves labelled 2 and 3 (4 and 5) in Figs. 4c, 5c and 6c (Figs. 4d, 5d and 6d)). 
The predicted spectra deformations in this case are rather small and, most probably, will 
be difficult to detect in SNO, Super Kamiokande and ICARUS experiments. Let us add 
that the distortions of the spectra can be enhanced by an appropriate choice of the specific 
normalization of the spectra, with the help of which one forms solar model independent spec- 
trum observables (as a comparison of Figs. 4 and 5 indicates). The results depicted in Figs. 
4 and 5 show that the vacuum oscillations, the MSW adiabatic, and the MSW nonadiabatic 
transitions of solar neutrinos lead to distinctly different deformations of the spectrum of the 
^B neutrinos, to be measured in the SNO (ICARUS) experiment. It seems very likely that 
the data from the SNO (ICARUS) detector on the ^B neutrino spectrum will allow one to 
test and to discriminate between these three possibilities. Adding the information about the 
seasonal variations of the signal will, most probably, permit to unambiguously distinguish 
between the vacuum oscillation and the MSW solutions. 

The distortions of the recoil-electron spectrum shown in Figs. 6 suggest that on the basis 
of the Super Kamiokande (ICARUS) data on the e^ spectrum alone it may be difficult to 
discriminate between vacuum oscillations with 4.5 x 10~^^eV^ < Am^ < 6.3 x 10~^^eV^ and 
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MSW nonadiabatic transitions of solar neutrinos (compare curves 5 and 6 in Fig. 6a, and 
2-5 in Fig. 6b with curves 4 in Fig. 6c and 1-3 in Fig. 6d). In the case of vacuum oscillations 
with Am^ from the above interval there will be seasonal variations of the signals in the Super 
Kamiokande, SNO and ICARUS detectors (see the dotted, dashed and long-dashed curves 
in Figs. 3d, 3e and 3f), which can be used to ehminate one of these two possibihties. 

It is also clear from Figs. 4, 5 and 6 that the information about the shapes of the *B 
neutrino and the e~ spectra to be obtained in the SNO, Super Kamiokande and ICARUS 
experiments, most probably, will not be sufficient to discriminate between an astrophysical 
and the large mixing angle (adiabatic) MSW solutions of the solar neutrino problem. How- 
ever, the measurement of the ratio of the number of events due to the solar neutrino induced 
charged current (CC) and neutral current (NC) reactions on deuterium, R*^^/^*^, to be per- 
formed with a relatively high precision in SNO experiment, will provide a crucial test of the 
large mixing angle MSW solution: for this solution one has R^^^^*^ = (0.3 — 0.4) Y{p'^/^'-' , 
where rCC'/^C' jg ^j^g value of the ratio predicted in the absence of oscillations and/or of 
MSW transitions. Note that the quantity Y{p^^^'^ does not depend on the total flux of 
neutrinos, and therefore is solar model independent; the value of Y{p^/^^ can be calculated 
with a high precision. 

To conclude, the envisaged capabilities of the next generation of solar neutrino experi- 
ments will allow one to perform crucial solar model independent tests of, and to discriminate 
between, the vacuum oscillation and the MSW solutions of the solar neutrino problem. It is 
very likely that the "solar neutrino puzzle" will be resolved by these experiments. 
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Figure Captions 

Fig. 1 The ratio of the vacuum oscillation {u^ t'^(r)) and of the SSM predicted signals 
(event rates) due to the ^B, pep, ^Be, and pp neutrinos as a function of the time of the year 
(in units of 365 days). The results shown are for experiments detecting the solar neutrinos 
via the u - e~ elastic scattering reaction (Super Kamiokande, BOREXINO, HELLAZ, etc.). 
The SSM predicted signals used represent the time independent one year average values of 
the event rates calculated within the model [5]. 

Fig. 2 The same as in Fig. 1 for different sets of values of Am^ and sin^ 29. Figures e and 
f correspond to i^s oscillations. 

Fig. 3 Time variations of the signals in Super Kamiokande (Figs. 3a, 3b, 3d and 3e) and 
SNO (Figs. 3d and 3f) detectors in the case of vacuum oscillations of ^B neutrinos, for 
the same values of the parameters Am^ and sin^26' for which Figs, la, lb,.... If and Figs. 
2a, 2b,..., 2f have been obtained (sohd, dotted, dashed, long-dashed, dash-dotted and long- 
dash-dotted lines in Figs. 3a-3c and 3d-3f, respectively). The signals are normahzed in such 
a way that in the absence of deviations from the standard 7% seasonal variation they will be 
constant in time and equal to 1 (i.e., horizontal lines crossing the vertical axis at the point 
1); the normalization used renders solar model independent the ratio of signals plotted. The 
results presented in Fig. 3a (3d) and Fig. 3b (3e) have been obtained with different values 
of the neutrino threshold energy: — 5 MeV and E^^ = 7.5 MeV, respectively. 

Fig. 4 Deformations of the '^B neutrino spectrum in the cases of i^g z/^(,-) (or z/g ^ Ug) 
oscillations (a and b), and of — J^ix{t) (or — > i^s) MSW transitions (c and d). The 
vacuum oscillation and the MSW spectra are divided by the SSM predicted spectrum [5], 
and each ratio of spectra is further normalized to the value this ratio has at E = 10 MeV. 
The double ratios plotted are solar model independent quantities. 

Fig. 5 The same as in Fig. 4, but with different normalization of the spectra (see eq. (18)). 
Fig. 6 Deformations of the spectrum of e~ from the reaction u + e~ — > u + e~ caused 



25 



by neutrinos, in the cases of oscillations in vacuum v^. <-> (a, and b (curves 1-4)), 
i^e ^ (b (curves 5 and 6)), and of MSW transitions v^, ^h(t) (c and d). Each of the 
predicted recoil-electron spectrum is divided by the standard one and the ratio so obtained 
is normalized to the value it has at Tg = 10 MeV. 

Fig. 7 Regions of values of the parameters Am^ and sin^26' allowed at 90% C.L. (dashed 
lines) and at 95% C.L. (sohd lines) by the current solar neutrino data in the case of MSW 
^fiir) (a and c) and Ug J^s (b and d) transitions of solar neutrinos. Figures a and 
b (c and d) have been obtained by including (without including) the uncertainties in the 
theoretical predictions [5] in the relevant x^— analysis. 
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